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A review of Big Bang Nucleosynthesis (BBN) is presented. Observations of deuterium and heIium-4 are dis- 
cussed. Some BBN restrictions on non-standard physics, especially on neutrino properties and time-variation of 
fundamental constants are given. 
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1. INTRODUCTION 

Big Bang Nucleosynthesis (BBN) is known to 
be one of three solid pillars on which the Standard 
Cosmological Model (SCM) stands. The other 
two include General Relativity (GR) and Cos- 
mic Microwave Background Radiation (CMBR). 
An agreement of BBN calculations of light ele- 
ment abundances with observations presents the 
strongest proof in favor of the statement that 12- 
14 billion years ago the universe was indeed hot 
with the temperatures in MeV range and that the 
entropy per baryon is huge, about 10^. 

According to the theory, light elements ^i/, 
^He, '^He, and ''Li have been created in the early 
universe during first few hundred seconds of her 
existence. The abundances of these elements span 
9 orders in magnitude and are in excellent, good, 
or reasonable agreement with the observational 
data, depending upon the moment when the com- 
parison of theory with the data is taken and upon 
the personal point of view of a researcher. 

The theory of BBN is robust, well defined, and 
quite precise. The largest uncertainty is intro- 
duced by the values of the cross-sections of nu- 
clear reactions. Theoretical accuracy is better 
than 0.1% for "Fe, better than 10% for '^H and 
is about 20-30% for Li ||l|,|. In all the cases 
theoretical uncertainty is much smaller than ob- 
servational precision. Observations of light ele- 
ments encounter two serious problems: system- 
atic errors and evolutionary effects. We will dis- 
cuss them below. 



In the next section physics of BBN and essen- 
tial parameters and inputs are described. In sec- 
tion 3 observational data are analyzed (looking 
from outside by a non-expert). In section 4 mod- 
ifications of the standard scenario are discussed. 
Conclusion is presented in the last section 5. 

2. PHYSICS OF BBN 

Physical processes which were essential for cos- 
mological creation of light elements took place at 
the temperatures, T, in the range of a few MeV 
down to 60-70 keV. The corresponding time, t, 
interval was from a fewxO.l sec up to 10'^ sec. At 
this time the universe was dominated by relativis- 
tic matter and the rate of cosmological cooling 
was determined by the expression: 
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where is the effective number of particle species 
in the cosmic plasma. In the standard model 
the factor 5, includes contributions from photons 
equal 2, from e^-pairs equal 7/2, and the contri- 
bution from three neutrino fiavors equal 3 • 7/4. 
Any additional, non-standard form of energy is 
parametrized in terms of effective number of neu- 
trino species: 

- 10.75 + J (A^. - 3) (2) 

Since the cosmological cooling rate depends upon 
it is clear that BBN is sensitive to any form 
of matter/ energy present in the primeval plasma 
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in the temperature interval (~ MeV) - 0.06 keV. 
This effect was noticed in refs. [||; detailed calcu- 
lations of the effect were pioneered in papers jj] . 

The first step in creation of light elements is 
"preparation" of neutrons. Their number density 
is determined by the reactions: 

n + ^ p + 

n + ^ p + e~ (3) 

Since the reaction rate is proportional to and 
the expansion rate is 77 ~ thermal equilibrium 
is maintained at high temperatures and n/p-ratio 
follows the equilibrium curve: 

n/p = e^p[~{dm + He) /T] (4) 

where Am — 1.293 MeV is the neutron-proton 
mass difference and /ie is the chemical potential 
of electronic neutrinos. In the standard model the 
latter is assumed to be vanishingly small. 

At T = 0.6 — 0.7 MeV expansion became faster 
than reactions ^ and the n/p-ratio would tend 
to a constant, if not the neutron decay. Be- 
cause of the decay the ratio slowly decreases as 
exp(— t/r„) with t„ = 889 sec. This behavior 
lasts approximately till T/vs = 65 keV when light 
element formation abruptly begins and all free 
neutrons quickly disappear from the plasma form- 
ing deuterium, helium, and lithium. The value 
of Tjvs is determined by the binding energies of 
light nuclei, in particular deuterium, and by the 
baryon-to-photon ratio, rj = nB/n^- 

Practically all neutrons that remained in the 
plasma to the moment when the temperature 
dropped down to Tjvs ended their lives in ^He 
because the latter has the largest binding energy. 
The mass fraction of the primordial ^iJe should 
be about 25%. A little deuterium and ^He sur- 
vived because they were not able to find a nu- 
cleon, N = p,n, for further transformation into 
^He through the chain of reactions: 

i^H, ^He) + N -^'^ He (5) 

As a result the abundance of primordial deu- 
terium would quickly decrease with rising baryon 
number density or, better to say, parameter 77, 
and its output with respect to hydrogen would 



be (a few) x 10~^ by number. A strong sensitivity 
of deuterium abundance to the baryon-to-photon 
ratio r] makes this element a very convenient tool 
to measure the baryonic charge of the universe. 
That's why primordial deuterium is called "bary- 
ometer" , the term suggested by David Schramm. 

An absence of a (quasi)stable nuclei with A = 5 
inhibits production of ^ii, because the latter 
should be produced by fusion of two lighter nu- 
clei and not in a sequence of reactions of free pro- 
tons or neutrons with nuclei. Correspondingly 
the fraction of '^Li is very small, (a few) x 10"^". 

The abundances of light elements as functions 
of the baryon number density, 7710 = lO^^ns/n^ 
are presented in fig. ^ taken from the review 

As we have seen, the abundances of light ele- 
ments essentially depend upon the following pa- 
rameters: 

1. Baryon-to-photon ratio, r]iQ — lO^'^UB/n^- 
Only two years ago this parameter was de- 
termined from BBN itself but now there is 
an independent way to measure it through 
the spectrum of angular fluctuations of 
CMBR. According to the BOOMERANG 
and DASI measurements: 

2_ r 0.021 ±0.004 [6], . 
"^'^ - 1 0.022tg:°- § , 

where fls = ps/pcrit is the fraction of 
baryon mass density in terms of cosmo- 
logical critical energy density and h = 
i?/100km/sec/Mpc is the dimensionless 
Hubble parameter, fis/i^ = 0.022 cor- 
responds to 7710 = 5.7. A higher value 
is given by MAXIMA-I measurements 
fis/i^ ^ 0.0325±0.007which disagrees with 
the other groups by two standard devia- 
tions. 

2. Rate of the reactions (H). It is expressed 
through neutron life-time which is pretty 
well known now, r„ = 886.7 ± 1.9 1^. 

3. Total cosmological energy density during 
BBN. This quantity is usually parametrized 
as the number of additional neutrino fam- 
ilies, AA^ = Ni, — 3. This is a precise 
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Figure 1. The predicted abundances of light ele- 
ments as a function of rjio ■ 



paranietrization if the additional energy is 
in the form of relativistic particles with the 
equation of state p — p/3. However it is 
not so for matter with a different equa- 
tion of state, e.g. nonrelativistic matter or 
vacuum-like energy. In this case the varia- 
tion of different abundances would be differ- 
ent from that induced by extra neutrinos. 

4. Neutrino degeneracy. It is assumed usually 
that the charge asymmetry of different neu- 
trino flavors is vanishingly small. Thus their 
chemical potentials are zero or negligible. 
Strictly speaking their values are unknown 
and the best way to determined them or ob- 



tain an upper bound on their magnitude is 
BBN. Chemical potentials of i^^ and i^r can 
be described by A7V, since their role is only 
to increase the energy density of relativis- 
tic matter at BBN. Element abundances are 
much more sensitive to chemical potential 
of Ve because the latter not only changes the 
canonical energy density but also directly 
shifts the rt/p-ratio (Q). 

3. OBSERVATIONAL DATA 

3.1. Observations of Deuterium 

The present day abundance of deuterium can 
be noticeably different from the primordial value 
because deuterium can be burnt in stars produc- 
ing ^He. Exact evolutionary effects are uncer- 
tain. One can rigorously say only that the ob- 
served today deuterium presents a lower bound 
to the primordial one. For more detail see e.g. re- 
view j5|. Fortunately in recent years observation 
of deuterium in large Z low metallicity (< 1/50 
of the solar) clouds of cold neutral gas (HI) be- 
came possible ^^Hl^. Such clouds are presum- 
ably not contaminated by stellar processes and 
the observed deuterium may be the primordial 
one; see however ref. |l5| where this statement is 
questioned. 

Since the deuterium isotope shift corresponds 
to velocity of only (-82) km/sec the clouds with a 
simple velocity structure are necessary for reliable 
data. In addition to the uncertainty induced by 
the unknown velocity structure, ionization correc- 
tions and a possible "interloper" further increase 
the systematic errors. A discussion of these ef- 
fects and a list of references can be found e.g. 
in |l|,0. 

Because of these problems, there are only a 
few measurements of deuterium at high red-shifts 
available now. They are summarized in Table |l|. 

The entries in Table |l| are situated in order of 
increasing the {D / H)-va\ues. Entries in columns 
3 and 5 correspond to the same system but the re- 
sult was changed to a larger value after its reanal- 
ysis in ref. ||2l[| where a complex velocity structure 
suggested by S II and H2 absorption lines was 
taken into account. 

The data seem to be grouped around three cen- 
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Table 1 

Abundance of deuterium in high red-shift Lyman alpha absorption systems 



10'^{D/H) - 1.5 



2.0762 |T7 


1 3.025 |9 


1 2.536 1 


2^ 


3.025 ^ 


1 3.572 1 


21 


2.504 H 


1 0.701 [ 


23 


1.65 ±0.35 




2.24 ±0.67 




2.54 ±0.23 


3.2 ±0.4 




3.25 ±0.3 


4.0 ±0.65 


20 ±5 





tral points. The majority of the data are accu- 
mulated around 10^ (D/H) = 3. There are two 
measurements indicating low deuterium about 1.5 
and one showing a huge value, 20 ± 5. Taken by 
the face value, there is an indication that there are 
systems with high, normal, and low deuterium 
abundances. If this is true, it would be a con- 
firmation of the model of ref. ||2^ where spa- 
tial variation of primordial abundances are pre- 
dicted. According to this work, two thirds of the 
sky should have normal abundance of deuterium, 
1/6 should have twice smaller ratio of D/H, and 
another 1 /6 should have about 5 times larger than 
normal abundance of deuterium. On the other 
hand, as is argued in ref. |26| , the line observed 
in the z = 0.701 absorption system |2^ which in- 
dicates a high deuterium abundance most proba- 
bly is not deuterium. This statement is based on 
a new measurement of the velocity dispersion of 
the neutral hydrogen. 

3.2. Observations of helium-4 

Helium-4 is very strongly binded nuclei so it 
cannot be destroyed by stellar processes but only 
produced in stars together with other heavier el- 
ements, "metals", e.g. oxygen (O) and nitrogen 
(N) - all elements heavier than *He are called 
"metals" in astronomy. Thus the observed mass 
fraction of '^He, should be larger than the pri- 
mordial one, Yp. ^He is observed in hot ionized 
gas regions (H II) in emission optical recombi- 
nation lines. In contrast to deuterium '^He has 
been observed only at relatively small distances, 
maximum, z — 0.045 |2^. This regions are con- 
taminated by stellar processes and to infer pri- 
mordial abundance from the data one should ex- 
trapolate to zero metallicity, i.e. to zero abun- 
dances of O or N. A detailed discussion of these 
issues can be found e.g. in the book |^^. The ex- 
trapolation to zero metallicity is realized by the 
linear function with the slope k = dY/dZ where 



Y is the observed mass fraction of '^He and Z is 
that of metals. Analysis by different groups give 
significantly different values n = 6.7 ± 2.3 ps| ], 
6.9 ± 1.5 1^, and k = 2.4 ± 1.0 [|l|. Though k 
is not accurately determined (or it may have dif- 
ferent values for different systems), the systems 
where '^He is studied have a low metallicity and 
the net effect on primordial abundance of *He is 
not very strong. 

The values of primordial mass fraction of '^He, 
Yp, measured by several different groups are pre- 
sented in table |[ The discrepancy between the 
results can be possibly prescribed to different 
treatment of correction factors: 1) ionization cor- 
rection factor (icf) which determines how much 
neutral hydrogen (invisible) is in the object under 
scrutiny, 2) temperature correction factor (tcf) 
which describes non-uniform temperature distri- 
bution, and 3) density of electrons. A recent 
analysis of the correction has been performed in 
ref. [Q, where it was shown in particular that 
combined icf and tcf would diminish the result 
by 0.002-0.004. It is worth noticing in this con- 
nection that the reanalysis of the data of ref. jsi]] 
in the paper where icf has been taken into ac- 
count, has led to a noticeably smaller result (see 
table I). 

Even the largest mass fraction of '^He pre- 
sented in table ^ corresponds to a rather low 
baryon number density, 7710 ~ 4.2 which is no- 
ticeably smaller than the value determined from 
CMBR (|). Possibly it means that systematic 
uncertainties and correction factors are underesti- 
mated and more work is necessary to obtain more 
precise results. This problem is discussed in the 
recent paper |^^. It seems that at the present 
time the accuracy in determination of primordial 
mass fraction of '^He is not sufhciently good. 
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Table 2 



Primordial mass fraction of ^i?e according to the results of different groups. 



Ref. 






30 


1 


32| 


@ 


1 


33 


1 




0.228 ±0.005 


0.234 ±0.002 


0.2348 ± 0.0025 


0.2443 ±0.0015 


0.238 ±0.003 



4. BIG BANG NUCLEOSYNTHESIS 
AND NON-STANDARD PHYSICS 

Many possible modifications of the standard 
cosmological scenario and/or Minimal Standard 
Model (MSM) in elementary particle physics can 
be strongly constrained or even excluded by BBN. 
As we have already mentioned, BBN permits to 
restrict the number of neutrino families. The 
most recent analysis |3^, based onYp — 0.238 ± 
0.002 ± 0.005, results in the following 95% CL 
upper limits: 

AiV < ( fo'- '^lo = 5.8 , 

\ 0.9 for ?7io = 2.4, ^ ' 

though in an earlier paper a much more re- 
strictive bound, AA'' < 0.2, was advocated. At 
the present stage is safe to assume that AA^ < 1. 
Hopefully in the near future one will be able to 
derive a stronger limit. 

Other additional parameters that can influence 
light element abundances are chemical potentials 
of different neutrino species, /ia, where a — e, fi, t. 
A possible role of neutrino degeneracy in big bang 
nucleosynthesis was noted already in the pioneer- 
ing work by Wagoner, Fowler, and Hoyle 
and after that it was analyzed in a number of 
papers. A combined analysis of the effect of si- 
multaneous variation of all three chemical poten- 
tials on BBN was recently performed in the pa- 
pers H^,^. Since the effect of positive ^e, that 
diminishes n/p-ratio, can be compensated by a 
non-zero /z^ or Ht, the limits obtained without 
any extra information are rather loose, assuming 
that the mentioned above conspiracy between /ie 
and iif^^r exists. However, additional data on the 
angular spectrum of CMBR permits to obtain 
stronger bounds even if the conspiracy is allowed: 

-O.OK^,, <0.2, |^,^,^,J <2.6 (8) 

These results are derived under assumptions that 
primordial fraction of deuterium is D/H = (3.0 ± 



0.4) • 10~^. Here £,a = fJ-a/T are dimensionless 
chemical potentials. 

These bounds are valid in absence of neutrino 
oscillations. On the other hand, solar and atmo- 
spheric neutrino anomalies present a strong evi- 
dence in favor of mixing between Vg, v^, and v-r 
(see e.g. talks by Y. Totsuka and A. McDonald 
at this Conference). Oscillations of these neutri- 
nos do not conserve individual leptonic charges 
and a primordial lepton asymmetries in elec- 
tronic, muonic, and tauonic charges would be re- 
distributed by the oscillations. In particular, a 
large /z^ or would create a noticeable fie of 
the same sign. An analysis of active-active neu- 
trino oscillations have not yet been accurately 
performed. Preliminary results |4l| are |^a| < 0.1 
for any flavor a — e, fj,,T. Let us note that in the 
case of vanishing chemical potentials oscillations 
between active neutrinos would not have a notice- 
able impact on BBN, because the neutrinos would 
practically remain in thermal equilibrium inde- 
pendently of oscillations and no deviation from 
the standard scenario would emerge. There could 
be a small effect due to spectral distortion of neu- 
trinos by late e'''e^-annihilation |^^. Oscillations 
could change the distorted spectrum and might 
be in principle observable in BBN. 

It is still not excluded that sterile companions 
of active neutrinos exist. Moreover, recent mea- 
surements [ p3[ seem to confirm creation of in 
beam. If this is the case, an interpretation of 
all neutrino data in terms of oscillations demands 
an existence of one (or several) sterile neutrino(s), 
Vs, mixed with active ones. 

In contrast to equilibrium active-active case, 
mixing of active and sterile neutrinos should have 
a noticeable impact on BBN. There are three pos- 
sible effects due to oscillations: 1) production of 
new relativistic states, i^s, and an increase of N^,; 
2) distortion of neutrino spectrum, especially im- 
portant is i^e spectrum; 3) resonance generation of 
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a large lepton asymmetry from a very small initial 
leptonic or baryonic ones (only the pioneer- 
ing and most recent papers are quoted; references 
to many other ones can be found therein). 

The effect of creation of additional relativis- 
tic species by oscillations is easily estimated and 
from the condition that BBN allows A7V extra 
neutrino families one can obtain the bound on 
vacuum mixing angle and the mass difference: 



(,5m VeV^) sin'' 261 < 



3.16- 10-5(AAf)2 
1.74- 10-^(A7V)2 



(9) 



for i/e Us and z^^.r Vg mixings respectively. This 
bounds are meaningful only if AA'^ < 1. 

Discussion of the effects of neutrino spectrum 
distortion on BBN can be found in rcf . |Q . The 
impact of asymmetry generation on BBN is dis- 
cussed in ref. [Q (second paper) and in rcfs. |46| . 
Resonance case is quite complicated and the ef- 
fect may have either sign. No simple bounds can 
be presented here. 

BBN permits put stringent bounds on possi- 
ble time- variation of fundamental constants. If 
we assume e.g. that the fine structure constant 
a is different at BBN from its present-day value, 
a = 1/137, we should expect that the neutron- 
proton mass difference should also change with 
time. The (n — p)-mass difference is given by 
dm = rrin — mp — nid — mu + aniem , where rud 
and m„ are the masses of u and d quarks and 
the last term describes electromagnetic loop con- 
tribution into Sm. As we have seen above the 
n/p-ratio is equal to n/p = exp(—Sm/Tf), where 
Tf is the freezing temperature of reactions (|^). 
This temperature is determined by the competi- 
tion between the Hubble expansion and weak in- 
teraction rates. The latter is proportional to the 
magnitude of the electroweak coupling constant, 
which is essentially a. Hence Tj ~ a^^^^. De- 
manding that successful results of BBN are not 
destroyed one can obtain that {Aa/a)BBN < (a 
few) X 10-2. 

5. Conclusion 

We see that gross features of BBN well agree 
with observations but the latter are not yet suffi- 
ciently accurate to make it really a precise science. 



Moreover there is a trend to discrepancy between 
the observations of deuterium which indicate a 
higher value of rjiQ than the observations of '^He. 
Hopefully it will be clear in a few coming years if 
this is a real problem or an artifact of systematic 
and evolutionary effects. We have not discussed 
above primordial Li because the accuracy of its 
measurements are rather low now but potentially 
this element could be very important for verifi- 
cation of the standard model. A recent discus- 
sion of '^Li can be found in ref. Still even 
with the existing level of accuracy BBN permit to 
put powerful constraints on deviations from the 
standard model. The number of extra neutrino 
species allowed by the contemporary observations 
is about unity and with this bound very little can 
said about mixing parameters between active and 
sterile neutrinos. However, if AA^ could be re- 
duced, say, to 0.3 the limit would be meaningful. 
The restriction on the time a is quite strong and 
may exclude some of the models predicting such 
variation. 
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